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t^! Abstract 
< 

£Nj ' A new functional calculus, developed recently for a fully non- 

perturbative treatment of quantum gravity, is used to begin a systematic 
construction of a quantum theory of geometry. Regulated operators corre- 
sponding to areas of 2-surfaces are introduced and shown to be self-adjoint 
\ on the underlying (kinematical) Hilbert space of states. It is shown that 

CN ' their spectra are purely discrete indicating that the underlying quantum 

geometry is far from what the continuum picture might suggest. Indeed, 
Q\ . the fundamental excitations of quantum geometry are 1-dimensional, rather 

like polymers, and the 3-dimensional continuum geometry emerges only on 
coarse graining. The full Hilbert space admits an orthonormal decompo- 
sition into finite dimensional sub-spaces which can be interpreted as the 
spaces of states of spin systems. Using this property, the complete spec- 
trum of the area operators is evaluated. The general framework constructed 
here will be used in a subsequent paper to discuss 3-dimensional geometric 



X 

operators, e.g., the ones corresponding to volumes of regions. 
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I. INTRODUCTION 



In his celebrated inaugural address, Riemann suggested |3| that geometry of space 
may be more than just a fiducial, mathematical entity serving as a passive stage for 
physical phenomena, and may in fact have a direct physical meaning in its own right. 
General relativity proved this vision to be correct: Einstein's equations put geometry on 
the same footing as matter. Now, the physics of this century has shown us that matter 
has constituents and the 3-dimensional objects we perceive as solids in fact have a discrete 
underlying structure. The continuum description of matter is an approximation which 
succeeds brilliantly in the macroscopic regime but fails hopelessly at the atomic scale. 
It is therefore natural to ask if the same is true of geometry. Does geometry also have 
constituents at the Planck scale? What are its atoms? Its elementary excitations? Is the 
space-time continuum only a "coarse-grained" approximation? If so, what is the nature 
of the underlying quantum geometry? 

To probe such issues, it is natural to look for hints in the procedures that have been 
successful in describing matter. Let us begin by asking what we mean by quantization of 
physical quantities. Let us take a simple example -the hydrogen atom. In this case, the 
answer is clear: while the basic observables -energy and angular momentum- take on a 
continuous range of values classically, in quantum mechanics their spectra are discrete. 
So, we can ask if the same is true of geometry. Classical geometrical observables such as 
areas of surfaces and volumes of regions can take on continuous values on the phase space 
of general relativity. Are the spectra of corresponding quantum operators discrete? If so, 
we would say that geometry is quantized. 

Thus, it is rather easy to pose the basic questions in a precise fashion. Indeed, they 
could have been formulated soon after the advent of quantum mechanics. Answering them, 
on the other hand, has proved to be surprisingly difficult. The main reason, it seems, is 
the inadequacy of the standard techniques. More precisely, the traditional approach to 
quantum field theory has been perturbative, where one begins with a continuum, back- 
ground geometry. It is then difficult to see how discreteness would arise in the spectra of 
geometric operators. To analyze such issues, one needs a fully non-perturbative approach: 
geometric operators have to be constructed ab-initio without assuming any background 
geometry. To probe the nature of quantum geometry, we can not begin by assuming the 
validity of the continuum picture. We must let quantum gravity itself decide whether 
this picture is adequate at the Planck scale; the theory itself should lead us to the correct 
microscopic picture of geometry. 

In this paper, we will use the non-perturbative, canonical approach to quantum gravity 
based on connections to probe these issues. Over the past three years, this approach has 
been put on a firm mathematical footing through the development of a new functional 
calculus on the space of gauge equivalent connections [4-11]. This calculus does not 
use any background fields (such as a metric) and is therefore well-suited to a fully non- 
perturbative treatment. The purpose of this paper is to use this framework to explore the 
nature of quantum geometry. 

In section 2, we recall the relevant results from the new functional calculus and out- 



2 



line the general strategy. In section 3, we present a regularization of the area operator. 
Its properties are discussed in section 4; in particular, we exhibit its entire spectrum. 
Our analysis is carried out in the "connection representation" and the discussion is self- 
contained. However, at a non-technical level, there is a close similarity between the basic 



ideas used here and those used in discussions based on the "loop representation" JT^JE 
Indeed, the development of the functional calculus which underlies this analysis itself was 
motivated, in a large measure, by the pioneering work on loop representation by Rovelli 



and Smolin [14]. The relation between various approaches will discussed in section 5. 



The main result of this paper should have ramifications on the statistical mechan- 



ical origin of the entropy of black holes along the lines of [T^.[T6|]. This issue is being 
investigated. 

II. PRELIMINARIES 

This section is divided into three parts. In the first, we will recall the basic 
structure of the quantum configuration space and, in the second, that of the Hilbert 



space of (kinematic) quantum states |T0| . The overall strategy will be summarized in the 
third part. 

A. Quantum configuration space 

In general relativity, one can regard the space AjQ of SU(2) connections modulo 
gauge transformations on a ("spatial") 3- manifold £ as the classical configuration space 



17H19[. For systems with only a finite number of degrees of freedom, the classical con- 
figuration space also serves as the domain space of quantum wave functions, i.e., as the 
quantum configuration space. For systems with an infinite number of degrees of free- 
dom, on the other hand, this is not true: generically, the quantum configuration space 
is an enlargement of the classical. In free field theory in Minkowski space (as well as 
exactly solvable models in low space-time dimensions), for example, while the classical 
configuration space can be built from suitably smooth fields, the quantum configuration 
space includes all (tempered) distributions. This is an important point because, typically, 
the classical configuration spaces are of zero measure; wave functions with support only 
on smooth configurations have zero norm! The overall situation is the same in general 
relativity. The quantum configuration space AjQ is a certain completion of AjQ ||||. 

The space AjQ inherits the quotient structure of AjQ, i.e., AjQ is the quotient of the 
space A of generalized connections by the space Q of generalized gauge transformations. 
To see the nature of the generalization involved, recall first that each smooth connection 
defines a holonomy along paths [Jin S: h p (A) := Pexp — J p A. Generalized connections 



1 For technical reasons, we will assume that all paths are analytic. An extension of the frame- 
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capture this notion. That is, each A in A can be defined [^,|J as a map which assigns 
to each oriented path p in E an element A(p) of SU(2) such that: i) A{p~ l ) = (A^p)) -1 ; 
and, ii) A(p 2 ° Pi) = A(p 2 ) ■ A(pi), where p~ l is obtained from p by simply reversing the 
orientation, p 2 o p 1 denotes the composition of the two paths (obtained by connecting 
the end of p\ with the beginning of p 2 ) and A(p 2 ) ■ A(pi) is the composition in SU(2). 
A generalized gauge transformation is a map g which assigns to each point v of E an 
SU(2) element g(x) (in an arbitrary, possibly discontinuous fashion). It acts on A in 
the expected manner, at the end points of paths: A{p) — > g>(t> + ) _1 ■ A(p) ■ g(v~), where 
V- and v + are respectively the beginning and the end point of p. If A happens to be a 
smooth connections, say A, we have A(p) = h p (A). However, in general, A(p) can not 
be expressed as a path ordered exponential of a smooth 1-form with values in the Lie 
algebra of SU(2) ||. Similarly, in general, a generalized gauge transformation can not be 
represented by a smooth group valued function on E. 

At first sight the spaces A, Q and A/Q seem too large to be mathematically con- 
trollable. However, they admit three characterizations which enables one to introduce 
differential and integral calculus on them We will conclude this sub-section by 

summarizing the characterization -as suitable limits of the corresponding spaces in lat- 
tice gauge theory- which will be most useful for the main body of this paper. 

We begin with some definitions. An edge is an oriented, 1-dimensional sub-manifold 
of E with two boundary points, called vertices, which is analytic everywhere, including 
the vertices. A graph in E is a collection of edges such that if two distinct edges meet, 
they do so only at vertices. In the physics terminology, one can think of a graph as a 
"floating lattice", i.e., a lattice whose edges are not required to be rectangular. (Indeed, 
they may even be non-trivially knotted!) Using the standard ideas from lattice gauge 
theory, we can construct the configuration space associated with the graph 7. Thus, we 
have the space Ay, each element A 7 of which assigns to every edge in 7 an element of 
SU(2) and the space Q 1 each element g 1 of which assigns to each vertex in 7 an element 
of SU(2). (Thus, if N is the number of edges in 7 and V the number of vertices, Ay 
is isomorphic with [SU(2)] N and Q 1 with [SU(2)] V ). Q 1 has the obvious action on Ay. 
Ary(e) — > g{v + )~ x ■ Ay(e) • giyJ). The (gauge invariant) configuration space associated 
with the floating lattice 7 is just A 7 /(? 7 . The spaces A, Q and A/Q can be obtained as 
well-defined (projective) limits of the spaces Ay, and Aj/Q 7 0,0- Note however that 
this limit is not the usual "continuum limit" of a lattice gauge theory in which one lets 
the edge length go to zero. Here, we are already in the continuum and have available to 
us all possible floating lattices from the beginning. We are just expressing the quantum 
configuration space of the continuum theory as a suitable limit of the configuration spaces 
of theories associated with all these lattices. 



work to allow for smooth paths is being carried out |2C|. The general expectation is that the 
main results will admit natural generalizations to the smooth category. In this article, A has the 
physical dimensions of a connection, (length) -1 and is thus related to the configuration variable 
A Q \d in the literature by A = GA Q \d where G is Newton's constant. 
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To summarize, the quantum configuration space A/Q is a specific extension of the 
classical configuration space A/Q. Quantum states can be expressed as complex- valued, 
square-integrable functions on A/Q, or, equivalently, as ^-invariant square-integrable 
functions on A. As in Minkowskian field theories, while A/Q is dense in A/Q topologically, 
measure theoretically it is generally sparse; typically, A/Q is contained in a subset set of 
zero measure of A/Q @. Consequently, what matters is the value of wave functions on 
"genuinely " generalized connections. In contrast with the usual Minkowskian situation, 
however, A, Q and A/Q are all compact spaces in their natural (Gel'fand) topologies [4-8]. 
This fact simplifies a number of technical issues. 

Our construction can be compared with the general framework of 'second quantization' 
proposed by Kijowski [^TJ] already twenty years ago. He introduced the space of states 
for a field theory by using the projective limit of spaces of states associated to a family 
of finite dimensional theories. He also, suggested, as an example, the lattice approach. 
The common element with the present approach is that in our case the space of measures 
on A is also the projective limit of the spaces of measures defined on finite dimensional 
spaces v4 7 . 



B. Hilbert space 



Since A/Q is compact, it admits regular (Borel, normalized) measures and for every 
such measure we can construct a Hilbert space of square-integrable functions. Thus, to 
construct the Hilbert space of quantum states, we need to select a specific measure on 
A/Q. 

It turns out that A admits a measure \i° that is preferred by both mathematical and 
physical considerations Mathematically, the measure \x° is natural because its defi- 

nition does not involve introduction of any additional structure: it is induced on A by the 
Haar measure on SU(2). More precisely, since A^ is isomorphic to [SU(2)] N , the Haar 
measure on SU(2) induces on it a measure fi° in the obvious fashion. As we vary 7, we 
obtain a family of measures which turn out to be compatible in an appropriate sense and 
therefore induce a measure fi° on A. This measure has the following attractive properties 
||: i) it is faithful; i.e., for any continuous, non-negative function / on A, J d\i° f > 0, 
equality holding if and only if / is identically zero; and, ii) it is invariant under the (in- 
duced) action of Diff[S], the diffeomorphism group of S. Finally, fi° induces a natural 
measure /2° on A/Q: pL° is simply the push-forward of \x° under the projection map that 
sends A to A/Q. Physically, the measure jx° is selected by the so-called "reality condi- 
tions". More precisely, the classical phase space admits an (over) complete set of naturally 
defined configuration and momentum variables which are real, and the requirement that 
the corresponding operators on the quantum Hilbert space be self-adjoint selects for us 
the measure p, |I(J . 



Thus, it is natural to use 7i° := L 2 (A/Q,dfi ) as our Hilbert space. Elements of 7i° 
are the kinematic states; we are yet to impose quantum constraints. Thus, 7i° is the 
classical analog of the full phase-space of quantum gravity (prior to the introduction of 
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the constraint sub- manifold) . Note that these quantum states can be regarded also as 
gauge invariant functions on A. In fact, since the spaces under consideration are compact 
and measures normalized, we can regard 7i° as the gauge invariant sub-space of the Hilbert 
space TC° := L 2 (A,dfi°) of square- integrable functions on A ||[7]]. In what follows, we we 
will often do so. 

What do "typical" quantum states look like? To provide an intuitive picture, we can 
proceed as follows. Fix a graph 7 with N edges and consider functions \I/ 7 of generalized 
connections of the form ^ 7 (A) = iJ)(A(ei), A(ejv)) for some smooth function ip on 
[SU(2)] N , where ei,...,ejv are the edges of the graph 7. Thus, the functions \I/ 7 know 
about what the generalized connections do only to those paths which constitute the edges 
of the graph 7; they are precisely the quantum states of the gauge theory associated 
with the "floating lattice" 7. This space of states, although infinite dimensional, is quite 
"small" in the sense that it corresponds to the Hilbert space associated with a system 
with only a finite number of degrees of freedom. However, if we vary 7 through all possible 
graphs, the collection of all states that results is very large. Indeed, one can show that it 
is dense in the Hilbert space TC°. (If we restrict ourselves to \l/ 7 which are gauge invariant, 
we obtain a dense sub-space in Ti°.) Since each of these states depends only on a finite 
number of variables, borrowing the terminology from the quantum theory of free fields 
in Minkowski space, they are called cylindrical functions and denoted by Cyl. Gauge 
invariant cylindrical functions represent the "typical" kinematic states. In many ways, 
Cyl is analogous to the space C^°(R 3 ) of smooth functions of compact support on R 3 
which is dense in the Hilbert space L 2 (R 3 , d 3 x) of quantum mechanics. Just as one often 
defines quantum operators - e.g., the position, the momentum and the Hamiltonians- on 
C£° first and then extends them to an appropriately larger domain in the Hilbert space 
L 2 (R 3 , d 3 x), we will define our operators first on Cyl and then extend them appropriately. 

Cylindrical functions provide considerable intuition about the nature of quantum 
states we are led to consider. These states represent 1-dimensional polymer-like exci- 
tations of geometry/gravity rather than 3-dimensional wavy undulations on flat space. 
Just as a polymer, although intrinsically 1-dimensional, exhibits 3-dimensional properties 
in sufficiently complex and densely packed configurations, the fundamental 1-dimensional 
excitations of geometry can be packed appropriately to provide a geometry which, when 
coarse-grained on scales much larger than the Planck length, lead us to continuum geome- 



tries [12], 22 1- Thus, in this description, gravitons can arise only as approximate notions 



in the low energy regime P3"| |. At the basic level, states in H° are fundamentally different 
from the Fock states of Minkowskian quantum field theories. The main reason is the 
underlying diffeomorphism invariance: In absence of a background geometry, it is not 
possible to introduce the familiar Gaussian measures and associated Fock spaces. 

C. Statement of the problem 

We can now outline the general strategy that will be followed in sections 4 and 5. 
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Recall that the classical configuration variable is an SU{2) connection^ A\ on a 6- 
manifold S, where % is the sw(2)-internal index with respect to a basis Tj. Its conjugate 
momentum E\ has the geometrical interpretation of an orthonormal triad with density 



weight one p4,17[) the precise Poisson brackets being: 



{Al(x),E<;(y)} = G5 b a S}5 3 (x,y), (2.1) 

where G is Newton's constant. (Recall from footnote 1 that the field A, used here, is 
related to A,id used in the literature via A = GA \d.) 

Therefore, geometrical observables -functionals of the 3-metric- can be expressed in 
terms of this field Ef. Fix within the 3-manifold £ any analytic, finite 2-surface S without 
boundary such that the closure of S in £ is a compact. The area As of S is a well-defined, 
real- valued function on the full phase space of general relativity (which happens to depend 
only on Ef). It is easy to verify that these kinematical observables can be expressed as: 

A s := I dx 1 A dx 2 \EfE 3i }^ , (2.2) 
Js 

where, for simplicity, we have used adapted coordinates such that S is given by x 3 = 0, 
and x 1 , x 2 parameterize S, and where the internal index % is raised by a the inner product 
we use on su(2), k(ri,Tj) = — 2Tr(7V7"j). 

Our task is to find the corresponding operators on the kinematical Hilbert space 7i° 
and investigate their properties. 

There are several factors that make this task difficult. Intuitively, one would expect 
that E°-(x) to be replaced by the "operator-valued distribution" — ihGS/SA l a (x). Unfortu- 
nately, the classical expression of As involves square-roots of products of E's and hence the 
formal expression of the corresponding operator is badly divergent. One must introduce 
a suitable regularization scheme. Unfortunately, we do not have at our disposal the usual 
machinery of Minkowskian field theories and even the precise rules that are to underlie 
such a regularization are not apriori clear. 

There are however certain basic expectations that we can use as guidelines: i) the 
resulting operators should be well-defined on a dense sub-space of H°; ii) their final 
expressions should be diffeomorphism covariant, and hence, in particular, independent of 
any background fields that may be used in the intermediate steps of the regularization 
procedure; and, iii) since the classical observables are real-valued, the operators should be 
self-adjoint. These expectations seem to be formidable at first. Indeed, these demands are 
rarely met even in Minkowskian field theories; in presence of interactions, it is extremely 
difficult to establish rigorously that physically interesting operators are well-defined and 



2 We assume that the underlying 3-manifold £ is orientable. Hence, principal SU(2) bundles 
over S are all topologically trivial. Therefore, we can represent the SU(2) connections on the 
bundle by a su(2)-valued 1-form on S. The matrices Tj are anti-Hermitian, given, e.g., by 
(— z/2)-times the Pauli matrices. 
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self-adjoint. As we will see, the reason why one can succeed in the present case is two- 
folds. First, the requirement of diffeomorphism covariance is a powerful restriction that 
severely limits the possibilities. Second, the background independent functional calculus 
is extremely well-suited for the problem and enables one to circumvent the various road 
blocks in subtle ways. 

Our general strategy will be following. We will define the regulated versions of area 
operators on the dense sub-space Cyl of cylindrical functions and show that they are 
essentially self-adjoint (i.e., admit unique self-adjoint extensions to This task is 

further simplified because the operators leave each sub-space H 7 spanned by cylindrical 
functions associated with any one graph 7 invariant. This in effect reduces the field theory 
problem (i.e., one with an infinite number of degrees of freedom) to a quantum mechanics 
problem (in which there are only a finite number of degrees of freedom). Finally, we 
will find that the operators in fact leave invariant certain finite dimensional sub-space 
of 7i° (associated with extended spin networks, introduced in Sec. 4.2). This powerful 
simplification further reduces the task of investigating the properties of these operators; 
in effect, the quantum mechanical problem (in which the Hilbert space is still infinite 
dimensional) is further simplified to a problem involving spin systems (where the Hilbert 
space is finite dimensional). It is because of these simplifications that a complete analysis 
is possible. 



III. REGULARIZATION 

Our task is to construct a well-defined operator A$ starting from the classical expres- 
sion ( |2.2| ). As is usual in quantum field theory, we will begin with the formal expression 
obtained by replacing Ef in ( |2.2j ) by the corresponding operator valued distribution Ef and 
then regulate it to obtain the required As- (For an early discussion of non-perturbative 
regularization, see, in particular, [p6| ). Our discussion will be divided in to two parts. 



In the first, we introduce the basic tools and, in the second, we apply them to obtain a 
well-defined operator As- 

To simplify the presentation, let us first assume that S is covered by a single chart 
of adapted coordinates. Extension to the general case is straightforward: one mimics the 
procedure used to define the integral of a differential form over a manifold. That is, one 
takes advantage of the coordinates invariance of the the resulting 'local' operator and uses 
a partition of unity. 



A. Tools 

The regularization procedure involves two main ingredients. We will begin by sum- 
marizing them. 

The first involves smearing of (the operator analog of) Ef(x) and point splitting of 



the integrand in ( |2.2j ). Since in this integrand, the point x lies on the 2-surface S, let 
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us try to use a 2-dimensional smearing function. Let f e (x,y) be a 1-parameter family of 
fields on S which tend to the S(x,y) as e tends to zero; i.e., such that 

lim/ d 2 y f e (x\x 2 ;y\y 2 )g(y\y 2 ) = g(x\x 2 ) , (3.1) 

for all smooth densities g of weight 1 and of compact support on S. (Thus, f e (x,y) is 
a density of weight 1 in i and a function in y.) The smeared version of Ef(x) will be 
defined to be: 

[Ef] f (x):= I d 2 yf e (x,y)Ef(y), (3.2) 
J s 

so that, as e tends to zero, [Ef]f tends to Ef(x). The point-splitting strategy now provides 
a 'regularized expression' of 

[A s ] f := ( d 2 x[[ d 2 yf e (x,y)Ef(y) [ d 2 z f € (x, z)E 3i (z)}? 
Js Js Js 

= [ d 2 x[[E!] f (x)[E%(x)]K (3.3) 
J s 

which will serve as the point of departure in the next subsection. To simplify technicalities, 
we will assume that the smearing field f e (x,y) has the following additional properties for 
sufficiently small e > 0: i) for any given y, f e (x, y) has compact support in x which shrinks 
uniformly to y; and, ii) f e (x,y) is non-negative. These conditions are very mild and we 
are thus left with a large class of regulators. 

We now introduce the second ingredient. To go over to the quantum theory, we want 
to replace Ef in ( |3.3| ) by Ef = —iGTi5/5A\. However, it is not apriori clear that, even after 
smearing, [Ef]f is a well-defined operator because: i)our wave functions \l/ are functionals 
of generalized connections A, whence it is not obvious what the functional derivative 
means; and, ii) we have smeared the operator only along two dimensions. Let us discuss 
these points one by one. 

First, let us fix a graph 7 and consider a cylindrical function \l/ 7 on A, 

%{A)=il>{A( ei ),..,A{e N )), (3.4) 

where, as before, N is the total number of edges of 7 and where ip is a smooth function 
on [SU(2)] N . Now, a key fact about generalized connections is that, for any given graph 
7, each A is equivalent to some smooth connection A 0: Given any A, there exists an A 
such that 

A(e k ) = h k [A] :=Pexp -/ A , (3.5) 



3 For example, f e (x,y) can be constructed as follows. Take any non-negative function / of 
compact support on S such that / d 2 xf(x) = 1 and set f e (x,y) = (l/e 2 )f((x — y)/e). Here, we 
have implicitly used the given chart to give f e (x,y) a density weight in x. 
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for all k = 1, N. (For any given A, the smooth connection A is of course not unique. 
However, this ambiguity does not affect the considerations that follow.) Hence, there is a 
1-1 correspondence between the cylindrical function \l/ 7 on A and function ip(hi(A), 
}ie(A)) on the space A of smooth connections and we can apply the operator [Ef]f to the 
latter. The result is: 

[E]] f (x) ■ %{A)=-iG%Y J J/yf^y){^^)y=o (^)(A) 

N ! 

= %e p l d 2 yf e (x,y)J2i dte I \t)5(y\e}(t))5(y^e 2 I (t))5(0,e 3 I (t)) 



1=1 



o 

dip 



(hj^tyhj^O^-^iA), (3.6) 
dh lB 

where, £p = V GTt is the Planck length, the index I labels the edges in the graph, [0, 1] 3 
t i — > e/(t) is any parameterization of an edge ej, hj(t',t) := V exp — f t A a (ej(s))e c }(s)ds 
is the holonomy of the connection A along the edge e/ from parameter value t tot'. Thus, 
the functional derivative has a well-defined action on cylindrical functions; the first of the 
two problems mentioned above has been overcome. 

However, because of the presence of the delta distributions, it is still not clear that 
[Ef]f is a genuine operator (rather than a distribution- valued operator). To explicitly see 
that it is, we need to specify some further details. Given a graph 7, we can just subdivide 
some of its its edges and thus obtain a graph 7' which occupies the same points in E as 
7 but has (trivially) more vertices and edges. Every function which is cylindrical with 
respect to the "smaller" graph 7 is obviously cylindrical with respect to the "larger" graph 
7' as well. The idea is to use this freedom to simplify the discussion by imposing some 
conditions on our graph 7. We will assume that: i) if an edge ej contains a segment which 
lies in S, then it lies entirely in the closure of S; ii) each isolated intersection of 7 with 
the 2-surface S is a vertex of 7; and, iii) each edge ej of 7 intersects S at most once. (The 
overlapping edges are often called edges 'tangential' to S; they should not be confused 
with edges which 'cross' S but whose tangent vector at the intersection point is tangent 
to S). If the given graph does not satisfy one or more of these conditions, we can obtain 
one which does simply by sub-dividing of some of the edges. Thus these conditions are 
not restrictive. They are introduced to simplify the "book-keeping" in calculations. 



Let us now return to ( p.6|) . If an edge ej has no point in common with S, it does 
not contribute to the sum. If it is contained in S, e| vanishes identically whence its 
contribution also vanishes. (For a subtlety, see the remark below Eq (|3.11| ).) We are 
thus left with edges which intersect S in isolated points. Let us first consider only those 
edges which are 'outgoing' at the intersection. Then, at the intersection point, the value 
of the parameter t is zero and, for a given edge ej, e'/ 3 is positive (negative) if ej is 
directed "upwards" along increasing a; 3 ("downwards" along decreasing x 3 ). Hence, ( |3.6| ) 
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becomes^: 



02 . N , ^ 



^ J 2 ^ 



,'/2 TV 

= E «/ fe& e,(0)) Li ■ V(A(ei), .... A( ejv )) , (3.7) 
A i=i 

where, the constant kj associated with the edge ej is given by: 

!0, if ej is tangential to S or does not intersect S, 
+ 1, if ej has an isolated intersection with S and lies above S (3.8) 
— 1, if ej has an isolated intersection with S and lies below S 

and where L\ is the left invariant vector field in the z-th internal direction on the copy of 
SU (2) corresponding to the 7-th edge 

L\ ■ iP(A( ei ), ...,A(e N )) = (A(ejy)i a , „ A . (3.9) 



d(A(ej)) 
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If some of the edges are 'incoming' at the intersection point, then the final expression of 
[Ef\f(x) can be written as: 

[Ef} f (x) ■ %= -^[^Kjf^v^Xj] ■ i>(A( ei ),...,A(e N )), (3.10) 
1 i=i 

where X\ is an operator assigned to a vertex v and an edge ej intersecting v by the 
following formula 

f(A(ei)T l ) 4 mz^^a , when e/ is outgoing 
-(^W)^, w hen e, is , n co mmg f" 1 ) 

Remark: Let us briefly return to the edges which are tangential to S. In this case, 
although Cj vanishes, we also have a singular term 5(0,0) (in the x 3 direction) in ( |3.6|) . 
Hence, to recover an unambiguous answer, for these edges, we need to smear also in 
the third direction using an additional regulator, say g t >(x 3 ,y 3 ). When this is done, one 
finds that the contribution of the tangential edges vanishes even before removing the 
regulator; as stated earlier, the tangential edges do not contribute. We did not introduce 
the smearing in the third direction right in the beginning to emphasize the point that this 
step is unnecessary for the edges whose contributions survive in the end. 

The right side again defines a cylindrical function based on the (same) graph 7. Denote 
by TC° the Hilbert space L 2 {A 1 ,dii°) of square integrable cylindrical functions associated 



4 In the first step, we have used the regularization dzg{z)5{z) = ^g{0) which follows if the 
5{z) is obtained, in the standard fashion, as a limit of functions which are symmetric about 0. 
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with a fixed graph 7. Since \i° is the induced Haar measure on A\ and since the operator 
is just a sum of right/left invariant vector fields, standard results in analysis imply that, 
with domain CyL of all C 1 cylindrical functions based on 7, it is an essentially self-adjoint 
on TC°. Now, it is straightforward to verify that the operators on 7i° obtained by varying 
7 are all compatible^ in the appropriate sense. Hence, it follows from the general results 
in m that [Ef]f(x), with domain Cyl 1 (the space of all C 1 cylindrical functions), is an 
essentially self-adjoint operator on TC°. For notational simplicity, we will denote its self- 
adjoint extension also by [Ef]f(x). (The context should make it clear whether we are 
referring to the essentially self-adjoint operator or its extension.) 

The fact that this operator is well-defined may seem surprising at first sight since 
we have used only a 2-dimensional smearing. Recall however that in free field theory 
in Minkowski space, the action of the momentum operator on cylindrical functions is 
well-defined in the same sense without any smearing at all. In our case, a 2-dimensional 
smearing is needed because our states contain one -rather than three- dimensional exci- 
tations. 

B. Area operators 

Let us now turn to the integrand of the smeared area operator (corresponding to ( |3.3| )). 
Denoting the determinant of the intrinsic metric on S by gs, we have: 

[g s ]f(x) ■ % := [E?] f (x)[E 3i ] f (x) ■ % 

= - e -f[T,< I ^)f^ v c ei )fe(x,v aj )X l I X)} ■ * 7J (3.12) 

4 I,J 

where the summation goes over all the oriented pairs (/, J); v ai and v aj are the vertices 
at which edges ej and ej intersect S; k(I, J) = KiKj equals if either of the two edges ej 
and ej fails to intersect S or lies entirely in S, +1 if they lie on the same side of S, and, — 1 
if they lie on the opposite sides. (For notational simplicity, from now on we shall not keep 
track of the position of the internal indices i\ as noted in Sec. 2.3, they are contracted 
using the invariant metric on the Lie algebra su(2).) The next step is to consider vertices 
v a at which 7 intersects S and simply rewrite the above sum by re-grouping terms by 
vertices. The result simplifies if we choose e sufficiently small so that, f e (x,v ai )f e (x,v aj ) 
is zero unless v„ T = v ni . We then have: 



5 Given two graphs, 7 and 7', we say that 7 > 7' if and only if every edge of 7' can be written as 
a composition of edges of 7. Given two such graphs, there is a projection map from to Ay, 
which, via pull-back, provides an unitary embedding of 7i°, into H°. A family of operators 
7 on the Hilbert spaces TC° is said to be compatible if U^^Oy = C? 7 Z7 7) y and U^^'Dy C D 7 
for all 7 > 7', where D 7 and Dy are the domains of 7 and O g > . 
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\gs] f (x) • * 7 = -\ [E(/ e (^««)) 2 E ^(4, W q X}J • * 7 , (3.13) 

a 11 

where the index a labels the vertices on S and I a and J a label the edges at the vertex a. 

The next step is to take the square-root of this expression. The same reasoning that 
established the self-adjoint ness of [Ef]f(x) now implies that [gs]f(x) is a non-negative self- 
adjoint operator and hence has a well-defined square-root which is also a positive definite 
self-adjoint operator. Since we have chosen e to be sufficiently small, for any given point 
x in S, f e (x,v a ) is non-zero for at most one vertex v a . We can therefore take the sum 
over a outside the square- root. One then obtains 

([&]/)*(*) • * 7 = f E Mx,v a )[ £ «(/ a , JjXlX'jJ* ■ * 7 - (3-14) 

a 11 

Note that the operator is neatly split; the x-dependence all resides in f € and the operator 
within the square-root is "internal" in the sense that it acts only on copies of SU(2). 

Finally, we can remove the regulator, i.e., take the limit as e tends to zero. By 
integrating both sides against test functions on S and then taking the limit, we conclude 
that the following equality holds in the distributional sense: 

H*) • *7 = % E * 2) (*,««)[ E <I a ,J*)X} a X^ ■ * 7 . (3.15) 

a 11 

Hence, the regularized area operator is given by: 

As ■ * 7 = f E [ E "Va, JJXiX}/* ■ * 7 . (3.16) 

(Here, as before, ct labels the vertices at which 7 intersects S and 7 Q labels the edges of 
7 at the vertex v a .) With Cyl 2 as its domain, As is essentially self-adjoint on the Hilbert 
space Ti°. 

Let us now remove the assumption that the surface £ is covered by a single chart of 
adapted coordinates. If such a global chart does not exist, we can cover £ with a family 
U of neighborhoods such that for each U G W there exists a local coordinates system (x a ) 
adapted to E. Let (<pu)ueu De a partition of unity associated to U. We just repeat the 
above regularization for a slightly modified classical surface area functional, namely for 



f dx 1 Adx 2 (pu[EfE 3i ]^ (3.17) 
J s 



which has support within a domain U of an adapted chart. Thus, we obtain the operator 
Asu- Then we just define 

A s = E Mu- (3.18) 
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The result is given again by the formula (|3.16 ). The reason why the functions ipu disappear 
from the result is that the operator obtained for a single domain of an adapted chart is 
insensitive on changes of this chart. This concludes our technical discussion. 

The classical expression As of ( [2.2|) is a rather complicated. It is therefore somewhat 
surprising that the corresponding quantum operators can be constructed rigorously and 
have quite manageable expressions. The essential reason is the underlying diffeomorphism 
invariance which severely restricts the possible operators. Given a surface and a graph, 
the only diffeomorphism invariant entities are the intersection vertices. Thus, a diffeo- 
morphism covariant operator can only involve structure at these vertices. In our case, it 
just acts on the copies of 577(2) associated with various edges at these vertices. 

We have presented this derivation in considerable detail to spell out all the assump- 
tions, to bring out the generality of the procedure and to illustrate how regularization can 
be carried out in a fully non-perturbative treatment. While one is free to introduce auxil- 
iary structures such as preferred charts or background fields in the intermediate steps, the 
final result must respect the underlying diffeomorphism invariance of the theory. These 
basic ideas will be used repeatedly for other geometric operators in the sequel to this 
paper. 



C. General properties of operators 

1. Discreteness of the spectrum: By inspection, it follows that the total area operator 
As leaves the sub-space of Cyl 2 which is associated with any one graph 7 invariant and 
is a self-adjoint operator on the sub-space 7i° of 7i° corresponding to 7. Next, recall that 
H" = L 2 (A~/, dn°), where A 7 is a compact manifold, isomorphic with (SU(2)) N where N 
is the total number of edges in 7. As, we explained below, the restriction of A s to H° 
is given by certain commuting elliptic differential operators on this compact manifold. 
Therefore, all its eigenvalues are discrete. Now suppose that the complete spectrum of 
As on 7i° has a continuous part. Denote by P c the associated projector. Then, given 
any \I/ in 7i°, P c • \P is orthogonal to 7i° for any graph 7, and hence to the space Cyl of 
cylindrical functions. Now, since Cyl 2 is dense in H.°, P c ■ \1/ must vanish for all \1/ in 7i°. 
Hence, the spectrum of A s has no continuous part. 

Note that this method is rather general: It can be used to show that any self-adjoint 
operator on H° which maps (the intersection of its domain with) 7i° to H°, and whose 
action on T~C° is given by elliptic differential operators, has a purely discrete spectrum on 
H.°. Geometrical operators, constructed purely from the triad field tend to satisfy these 
properties. 

2. Area element: Note that not only is the total area operator well-defined, but in 
fact it arises from a local area element, which is an operator- valued distribution in 
the usual sense. Thus, if we integrate it against test functions, the operator is densely 
defined on 7i° (with C 2 cylindrical functions as domain) and the matrix elements 

(*V> jt s (x) ■ %) (3.19) 
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are 2-dimensional distributions on S. Furthermore, since we did not have to renormalize 
the regularized operator ( p.!4j ) before removing the regulator, there are no free renormal- 
ization constants involved. The local operator is completely unambiguous. 

3. [gs\f versus its square-root: Although the regulated operator [g s ]f is well-defined, 
if we let e to go zero, the resulting operator is in fact divergent: roughly, it would lead 
to the square of the 2-dimensional 5 distribution. Thus, the determinant of the 2-metric 
is not a well-defined in the quantum theory. As we saw, however, the square-root of the 
determinant is well defined: We have to first take the square-root of the regulated expres- 
sion and then remove the regulator. This, in effect, is the essence of the regularization 
procedure. 

To get around this divergence of gs, as is common in Minkowskian field theories, we 
could have first rescaled [g,s}f] by an appropriate factor and then taken the limit. Then 
result can be a well-defined operator, but it will depend on the choice of the regulator, 
i.e., additional structure introduced in the procedure. Indeed, if the resulting operator 
is to have the same density character as its classical analog gs{x) -which is a scalar 
density of weight two- then the operator can not respect the underlying diffeomorphism 
invariance.0 For, there is no metric/chart independent distribution on 5* of density weight 
two. Hence, such a 'renormalized' operator is not useful to a fully non-perturbative 
approach. For the square-root, on the other hand, we need a local density of weight one. 
And, the 2-dimensional Dirac distribution provides this; now is no apriori obstruction for a 
satisfactory operator corresponding to the area element to exist. This is an illustration of 
what appears to be typical in non-perturbative approaches to quantum gravity: Either the 
limit of the operator exists as the regulator is removed without the need of renormalization 
or it inherits back-ground dependent renormalization fields (rather than constants). 

4. Vertex operators: As noted already, in the final expressions of the area element and 
area operators, there is a clean separation between the 'x-dependent' and the 'internal' 
parts. Given a graph 7, the internal part is a sum of square-roots of the operators 

As tBb := E <I«, JJXj.Xj. (3.20) 

la ■ J a. 

associated with the surface S and the vertex v a on it. It is straightforward to check 
that operators corresponding to different vertices commute. Therefore, to analyze the 
properties of area operators, we can focus just on one vertex operator at a time. 

Furthermore, given the surface S and a point v on it, we can define an operator A s, v 
on the dense sub-space Cyl 2 on H° as follows: 



6 If, on the other hand, for some reason, we are willing to allow the limiting operator to have a 
different density character than its classical analog, one can renormalize [g]f(x) in such a way as 
to obtain a background independent limit. For instance, we may use f e = (l/e 2 )9(\x — x'\ — |), 
and rescale [g]f by e 2 before taking the limit. Then the limit is a well defined, diffeomorphism 
covariant operator but it is a scalar density of weight one rather than two. 
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A xv — / K ( J ' J ) X i X J ■ if 7 intersects 5 in u, , . 

A ^'^-"lO, Otherwise 

where 7 and J label the edges of 7 which have v as a vertex. (Recall that every cylindrical 
function is associated with some graph 7. As before, if 7 intersects S at v but v is not 
a vertex of 7, one can extend 7 just by adding a new vertex v and orienting the edges 
at v to outgoing.) It is straightforward to verify that this definition is unambiguous: if 
a cylindrical function can be represented in two different ways, say as \P 7 and \I/y, then 
&s,v " ^7 and A S)V ■ \l/y are two representations of the same function on A. There is a 
precise sense in which As >v can be regarded as a Laplacian operator on TC°. The area 
operator is a sum over all the points v of S of square-roots of Laplacians, 



(Here the sum is well defined because, for any cylindrical function, it contains only a 
finite number of non-zero terms, corresponding to the isolated intersection points of the 
associated graph with S). We will see in the next subsection that this fact is reflected in 
its spectrum. 

5. Gauge invariance: The classical area element ^J~g~s is invariant under the internal 
rotations of triads Ef; its Poisson bracket with the Gauss constraint functional vanishes. 
This symmetry is preserved in the quantum theory: the quantum operator ^fg~s commutes 
with the induced action of Q on the Hilbert space 7i° . Thus, y/g~s and the total area 
operator A s map the space of gauge invariant states to itself; they project down to the 
Hilbert space 7i° of kinematic states. 

Note, however, that the regulated triad operators [Ef]f are not gauge invariant; they 
are defined only on 7i°. Nonetheless, they are useful; they feature in an important way 
in our regularization scheme. In the loop representation, by contrast, one can only in- 
troduce gauge invariant operators and hence the regulated triad operators do not exist. 
Furthermore, even in the definition (|3.3|) of the regularized area element, one must use 
holonomies to transport indices between the two points y and z. While this manifest 
gauge invariance is conceptually pleasing, in practice it often makes the calculations in 
the loop representation cumbersome; one has to keep track of these holonomy insertions 
in the intermediate steps although they do not contribute to the final result. 

6. Overall Factors: The overall numerical factors in the expressions of various op- 
erators considered above depend on two conventions. The first is the convention noted 
in footnote 4 used in the regularization procedure. Could we not have used a different 
convention, setting / °° dzg(z)5(z) = cg(0) and J oo dzg(z)5(z) = (1 — c)g(0) for some 
constant c 7^ 1/2? The answer is in the negative. For, in this case, the constant Ki would 
take values: 

0, if ej is tangential to S or does not intersect S, 

ki = \ +2c, if ei has an isolated intersection with S and lies above S (3.23) 

-2(1 — c), if e/ has an isolated intersection with S and lies below S 
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It then follows that, unless c = 1/2, the action of the area operator A$ on a given 
cylindrical function would change if we simply reverse the orientation on S (keeping the 
orientation on E the same). Since this is physically inadmissible, we must have c = 1/2; 
there is really no freedom in this part of the regularization procedure. 

The second convention has to do with the overall numerical factor in the action, which 
dictates the numerical coefficients in the symplectic structure. Here, we have adopted 
the convention of (see chapter 9) which makes the Poisson bracket {A % a (x), Ej(y)} = 
G5l5 l j5(x, y), enabling us to express Ef(x) as —iGTiS / 5 A % a (x) . (Had we rescaled the action 
by 1/87T as is sometimes done, in our expressions, Newton's constant G would be replaced 
by 87rG.) 

IV. EIGENVALUES AND EIGENVECTORS 

This section is divided into three parts. In the first, we derive the complete spectrum 
of the area operators, in the second, we extend the notion of spin networks, and in the 
third, we use this extension to discuss eigenvectors. 

A. The complete spectrum 

We are now ready to calculate the complete spectrum of As- Since As is a sum of 
square-roots of vertex operators which all commute with one another, the task reduces 
to that of finding the spectrum of each vertex operator. Furthermore, since vertex op- 
erators map (C 2 ) cylindrical functions associated with any one graph to (C°) cylindrical 
functions associated with the same graph, we can begin with an arbitrary but fixed graph 
7. Consider then a vertex operator A 5^ and focus on the edges of 7 which intersect S 
at v. Let us divide the edges in to three categories: let ei,...,e<j lie 'below' S ('down'), 
ea+i, e u lie 'above' S ('up') and let e u+ i, ...,e t be tangential to S. (As before, the labels 
'down' and 'up' do not have an invariant significance; the orientation of S and of E enable 
us to divide the non-tangential edges in to two parts and we just label one as 'down' and 
the other as 'up'.) Let us set: 

jfl 7 := ~* + ••• + ^Q)> <4> 7 := ~i PQ+i + ••• + -^u)> 

j{t) 1 ,-/yt I , yi\ j(d+u) 1 j(d) 1 j(u) 1 / 41 x 

J S,v 7 ~ % l A «+l + ••• + A J' J S,v J S,v + J S,v l 4 - 1 ^ 

where X) is the operator defined in ( pMl| ) assigned to the point v and an edge ej at v. 
This notation is suggestive. We can associate with each edge e a particle with only a spin 
degree of freedom. Then, the operators — iX l e can be thought of as the ith component 

of angular momentum operators associated with that particle and J s d \ , jj^l and J$\ as 
the total 'down', 'up' and 'tangential' angular momentum operators at the vertex v. 
By varying the graph, we thus obtain a family of operators. It is easy to check that 

they satisfy the compatibility conditions and thus define operators J s u \ , J s d \ .Jg\ and 
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jjtv^ on Cyl. It is also easy to verify that they all commute with one another. Hence, 
one can express the vertex operator As )t) simply as: 

-As,, = (45* -45') (45* -45'); 

because of the factor k(I, J) in ( |3.21| ), the edges which are tangential do not feature in 
this expression. 

The evaluation of possible eigenvalues is now straightforward. It is simplest to express 

A SjV as 

- As,, = 2(4J) 2 + 2( Jg) 2 - (jiT') 2 ■ (4-3) 

and^ as in elementary text-books, go to the representation in which the operators 
(4>) 2 > (4>) 2 an d (4«" ) 2 are diagonal. If we restrict now the operators to Cyl 7 as- 
sociated to a fixed graph, it is obvious that the possible eigenvalues A of As, v are given 
by: 

X s , v = 2jW(j® + 1) + 2jM(j«> + 1) - + 1) (4.4) 

where j^,j( u ' and j( d+n ) are half integers subject to the usual condition: 

3 {d+u) e {\j (d) -J {u) \, \j (d) -j (u) \ + 1, - ,J {d) +J (u) }- (4.5) 

Returning to the total area operator, we note that the vertex operators associated 
with distinct vertices commute. Although the sum ( 3.22|) is not finite, restricted to any 
graph 7 and Cyl 7 it becomes finite. Therefore, the eigenvalues as of As are given by 

*s = % E H d) Ui d) + 1) + V£W + 1) - 3i d+u) U { « +u) + 1)] ^ (4-6) 

where a labels a finite set of points in S and the non-negative half-integers assigned to 
each a are subject to the inequality (|4.5|) . The question now is if all these eigenvalues are 
actually attained, i.e., if, given any as of the form (|4.6|), there are eigenvectors in H° with 
that eigenvalue. In Sec. 4.3, will show that the full spectrum is indeed realized on 7i°. 

The area operators map the subspace H° of gauge invariant elements of 7i° to itself. 
Hence we can ask for their spectrum on 7i°. We will see in Sec. 4.3 that further restrictions 
can now arise depending on the topology of the surface S. There are three cases: 

• i) The case when S is an open surface whose closure is contained in E. An example 
is provided by the disk z = 0, x 2 + y 2 < r in R 3 . In this case, there is no additional 
condition; all a s of fl4.6p subject to fl4.5[) are realized. 



ii) The case when the surface S is closed (dS = 0) and divides E into disjoint open 
sets Ei and E 2 (i.e., E = Ei U S U E 2 with Ei fl S 2 = 0.) An example is given by: 
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X = R 3 and S = S 2 . In this case, there is an a condition on the half integers 
and j^f 1 that appear in (|4.6[) in addition to 



£jl d) =iV, an d ]TjW = iV' ( 4 -7) 



for some integers N and N' . 

• hi) The case when 5* is closed but not of type ii). An example is given by; X = 
S 1 x S 1 x S 1 and S = S 1 x S 1 . In this case, the additional condition is milder: 

EVP + = N ( 4 - 8 ) 

a 

for some interger N. 

Next, let us note some properties of this spectrum of Ag. By inspection, it is clear 
that the smallest eigenvalue is and that the spectrum is unbounded from above. One 
can ask for the 'area gap' i.e., the value of the smallest non-zero eigenvalue. On the full 
Hilbert space H°, it is given by: 

a° s = ^4- (4.9) 



This is a special case of the situation when there is only one term in the sum in ( |4.6|) with 

j( d ) -0 jW= j(d+u) _ j_ Then 



as = f \fj(j + l), (4-10) 

and, if we choose j — |, we obtain the eigenvalue a° s . On the Hilbert space 7i° of gauge 
invariant states, on the other hand, because of the constraints on the spectrum discussed 
above, the area gap is sensitive to the topology of S: 

a° s = —T-£p if S is of type i) 
a° s = —^-£ P if S is of type ii) 

a° s = -£% if S is of type hi). (4.11) 

Another important feature of the spectrum is its behavior for large as- As noted above, 
the spectrum is discrete. However, an interesting question is if it approaches continuum 
and, if so, in what manner. We will now show that as as —>■ oo, the difference Aas 
between as and its closest eigenvalue satisfies the inequality 

Aa s < (4/2)(VV^) + 0((4/as))4 (4.12) 
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and hence tends to zero (irrespective of the topology of S). Specifically, given (odd) 
integers M and N satisfying 1 < M < 2^/N, we will obtain an eigenvalue as,N,M of As 
such that for sufficiently large N, the bound (|4.12|) is explicitly realized. [] Let us label 



representations of SU(2) by their dimension, n a = 2j a + 1. Let n a , a = 1,...,M be 
(odd) integers such that Z) a=1 n a = N, and \n a — jj\ < 2 Then, for each M, we have 
from (|4.10|) an eigenvalue as. 



N,M 



1)1 M 
-P 



a S ,N,M = -7T^2\jia{ia + l) 
Z a=l 



02 M 



V-£ + *£ + 0(i)) (4.13) 



4 v 2N N 2 X N' 

for some integer k E [1, M/2]. As M varies between 1 and 2^/N, as t N,M varies between 
(£ 2 P /A)N and (£ 2 P /A)(N-2)+4k/N < (£ 2 P /4)((N -2)+4/^N). Hence', given a sufficiently 
large as, there exist integers N, M satisfying the conditions given above such that Aas '■ = 
I as — a n,m\ satisfies the inequality (|4.12|) . 

We will conclude this discussion of the spectrum by providing an alternative form of 
the expression ( |4.4|) which holds for gauge invariant states. This form will be useful in 
comparing our result with those obtained in the loop representation (where, from the 
beginning, one restricts oneself to gauge invariant states.) Let \l/ 7 be a gauge invariant 
cylindrical function on A. Then, the Gauss constraint implies that, at every vertex v of 
7, the following condition must hold: 

Y,Xj-% = Q, (4.14) 

where / labels the edges of 7 at the vertex v and X\ is assigned to the point v and vertex 
ej (see ( |3.11| ). Therefore, 



48 + 45 + Jgi" = 0. (4.i5) 

Hence, one can now express the operator ( |4.3| ) in an alternate form, 

-A s , v = 2(J s y + 2(J s y-(J s y. (4.16) 



7 This calculation was motivated by the results of Bekenstein and Mukhanov |l5| and our 
estimate has an interesting implication on whether the Hawking spectrum is significantly altered 
due to quantum gravity effects. Because the "level spacing" Aas goes to zero as as goes to 
infinity, the considerations of jyj do not apply to large black holes in our approach and there 
is no reason to expect deviations from Hawking's semi-classical results. On the other hand, for 
small black-holes, -i.e., the final stages of evaporation- the estimate does not apply and one 
expects transitions between area eigenstates to show significant deviations. 
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Furthermore, if it happens that 7 has no edges which are tangential to S at v, ( 4.14j) 
implies: 



As,, = 4(J^) 2 = 4(4!) 2 , (4.17) 



whence the corresponding restricted eigenvalues of A$ are given by J2^pyj(j + 1)> where 
j are half-integers. 



B. Extended spin networks 

As a prelude to the discussion on eigenvectors, in this sub-section we will generalize 
the constructions and results obtained in |],|1(],[?7]] on spin networks and spin network 
states. The previous work showed that the spin network states provide us with a natural 
orthogonal decomposition of the Hilbert space 7i° of gauge invariant states in to finite 
dimensional sub-spaces. Here, we will extend those results to the space H°. 

We begin by fixing some terminology. Given N irreducible representations 71*1, 7Tjv of 
577(2), an associated invariant tensor c m k+i--- m N mi mk is a multi-linear map from (g>/=i 717 
to (g)f =k+1 Uj such that: 

WS)^---^^ = c"*— ni ... nk , (4.18) 

for arbitrary g G SU(2), where 717(g) is the matrix representing g in the representation 
717. An invariant tensor c mi ,„ mk mk+1 '" mN is also called an intertwining tensor from the 
representations 717, ...,7Tfc into Tr^+i, tt^. All the invariant tensors are given by the 
standard Clebsch- Gordon theory. 

An extended spin network is an quintuplet (7, 7?, c, p, M) consisting of: 

i) A graph 7; 

ii) A labeling 7? := (717, .., 7i7v) of the edges ei,...,ejv of that graph 7 with irreducible 
and non-trivial representations of SU(2); 

iii) A labeling p : = (pi, ...,/v) of the vertices vi, ...,ty of 7 with irreducible representa- 
tions of SU(2), the constraint being that for every vertex v a the representation p a 
emerges in the decomposition of the tensor product of representations assigned by 
7T to the edges intersecting v a ; 

iv) A labeling c = (ci, cy) of the vertices Vi, vy of 7 with certain invariant tensors, 
namely, assigned to a vertex v a is an intertwining tensor c a from the representations 
assigned to the edges coming to v a and p a to the representations assigned to the 
outgoing edges at v a ; and, 

v) A labeling M := (M a ) a= i t _ t v = (Mi, M v ) of the vertices Vi,...,Vy of 7 which 
assigns to every vertex v a a vector M a in the representation p a \ 
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It should be emphasized that every ttj is necessarily non-trivial whereas p a may be 
trivial (i.e., 1-dimensional). In the gauge invariant context [|9,10], p a are all trivial whence 



the items iii) and v) are unnecessary. The details of these conditions may seem somewhat 



complicated but they are necessary to achieve the orthogonal decomposition (|4.22j ). 

From spin networks, we can construct states in Ti,°. An extended spin network state 
J\f y £M is simply a C°° cylindrical function on A constructed from an extended spin net- 
work (7, 7?, p, c, M), 

N V 

KcXi W ■= [® MM*)) ® <g) M a ] ■ [ ®Li c a ], (4.19) 

1=1 a=l 

for all A £ A, where, as before, A(e{) is an element of G associated with an edge ej 
and '•' stands for contracting, at each vertex v a of 7, the upper indices of the matrices 
corresponding to all the incoming edges, the lower indices of the matrices assigned to 
all the outgoing edges and the upper index of the vector M a with all the corresponding 
indices of c a . (We skip 7? and p in the symbol for the extended spin network function 
because the intertwiners c contain this information.) Thus, for example, in the simple 
case when the network has only two vertices, and all edges originate at the first vertex 
and end at the second, A/" 7?A j can be written out explicitly as: 

KcM = MMei))Z -. MM**))Z Mf'M?* cT- m » mll c 2ni ... nNm>2 , (4.20) 

where indices mi, nj range over 1, 2ji + 1 and m' a ranges over 1, 2j a+ \. Given any 
spin network, ( [4.19j ) provides a function on A which is square-integrable with respect to 
the measure pf . Given an extended spin network function on A, the the range R(j) of 
the associated graph 7 is completely determined. Thus, two spin networks can define 
the same function on A if one can be obtained from the other by subdividing edges and 
changing arbitrarily the orientations. 

It turns out that the spin network states provide a decomposition of the full Hilbert 
space TC° into finite dimensional orthogonal sub-spaces (compare with PJlQfl). Given a 
triplet (7, 7r, p) defined by (i — iii) above, consider the vector space Ti. 7 ^^ spanned by the 
spin network functions A/V ,7 given by all the possible choices for c, M compatible with 
fixed labelings tt, p. Note that, according to the representation theory of compact groups, 
every H. 1 ^ % p is a finite dimensional irreducible representation of Q in Cyl. The group acts 
there via 

KcM^Ag) = ^^(A), M' a = Pa (g(v a ))M a . (4.21) 
Modulo the obvious completions, we have the following orthogonal decomposition, 

^° = Htm (4.22) 

where, given a graph 7, the labelings 7r and p range over all the data defined above by 
(i-iii) whereas for 7 in the sum we take exactly one representative from every range of an 
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analytic graph in E. When p is trivial we skip p in Ji 1 ^ i p. On 7i 7j #, the action of the gauge 
transformations group Q is trivial and we have the following orthogonal decomposition of 
the Hilbert space of gauge invariant cylindrical functions, 

^° = Kr,*, (4-23) 

fl( 7 ),7f 



where we used the same conventions as in ([4.22 ) . Thus, we recover the result on spin 
network states obtained in PJlO]. 

We conclude this sub-section with a general comment on spin network states. Consider 
tri-valent graphs, i.e. graphs 7 each vertex of which has three (or less) edges. In this case, 
the standard Clebsch- Gordon theory implies that the the number of associated gauge 
invariant spin network functions is severely limited: the corresponding sub-space of 7i° 
is one dimensional. Hence, on the sub-space Cyl of 7i° corresponding only to tri-valent 
graphs, the (normalized) spin network states provide a natural orthonormal basis. What is 
remarkable is that these spin networks were first introduced by Penrose f28[ already twenty 



five years ago to probe the microscopic structure of geometry, although in a different 
context. Because of the simplicity (and other attractive properties) of these Penrose spin 
network states it is tempting to hope that they might suffice also in the present approach 
to quantum gravity. Indeed, there were conjectures that the higher valent graphs are 
physically redundant. However, it turns out that detailed physical considerations rule out 
this possibility; quantum gravity seems to need graphs with unlimited complexity. 



C. Eigenvectors 

We are now ready to exhibit eigenvectors of the operators As„ and A$ for any of 
the potential eigenvalues found in section 4.1. We will begin with the full, non-gauge 
invariant Hilbert space 7i° and consider an arbitrary surface S. Since 7i° serves as the 
(gravitational part) of the kinematical Hilbert space in theories in which gravity is coupled 
to spinor fields, our construction is relevant to that case. In the second part of this sub- 
section, we will turn to the gauge invariant Hilbert space 7i° and exhibit eigenvectors for 
the restricted range of eigenvalues presented in section IV A . 



Fix a point v in the surface S. We will investigate the action of the operators (jjfl) 2 , 

{J^l) 2 , {Jgv Y and As, v on extended spin network states. Without loss of generality 
we can restrict ourselves to graphs which are adapted to S and contain v as a vertex, say 
v = v\. Given a graph 7 and labeling 7? and p of its edges and vertices by representations 
of SU(2), we shall denote by C v the linear space of the intertwining tensors which are 



compatible with 7? and p at v in the sense of section [IV B| . Let (7,7?, p,c,M) be an 
extended spin network and J\l ' SJ ^ be the corresponding state. As one can see from Eqs. 



|4.1| , |3.21|) , each of the four operators above is given by a linear combination (with constant 



coefficients) of gauge invariant terms of the form bi 1 ___i E XY 1 ---X\ E where 6^...^ is a constant 
tensor and all the Xs are associated with the point v and the edges which meet there. 
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On A/V - r, the action of any operator of this type reduces to a linear operator o v acting 
in C v . More precisely, if O is any of the above operators, we have 

0^*0 = Mr** ( 424 ) 

where A/" 7?A j is again an extended spin network state and the network (7, jf, p, <?, M) 
differs from the first one only in one entry of the labeling <? corresponding to the vertex 
v; c' a = c a for all the vertices v a 7^ v and c\ = o v c\. Consequently, the problem of 
diagonalizing these operators reduces to that of diagonalizing a finite symmetric matrix 
of o v . Note that a constant vector M assigned to v does not play any role in this action 
and hence will just make eigenvectors degenerate. 

In the case of operators (jjfl) 2 , (jjsl) 2 and (Jg 1 ^) 2 , the (simultaneous) eigenstates are 
given by the group representation theory. We can now spell out the general construction. 

Let us fix a graph 7 and arrange the edges that meet at v into three classes as be- 
fore; ei,...,ed] Cd+i, e u ; e u+i ,...,e t . Let us also fix a labeling iri,...,ir t of these edges 
by irreducible, non-trivial representations of SU(2) and an irreducible (possibly trivial) 
representation p which emerges in the decomposition of 7Ti £5 ... £5 iTf Consider now the 
following ingredients: 

a) Irreducible representations p(d), P{u) and p(d+ u y, 

b) Invariant tensors c^ mi - mdm ' , c ™rf+i- m " m anc [ C( u+d ) m 'm" m associated, respec- 
tively, to the representations 7Ti, 71^, //(<*), and to n d+1 , ...,n u , p( u ) and finally to 
P(d), P(u), P(d+u)', and, 

c) Invariant tensor c ^ n mu+1 "- mtrn assoc i a ted to p( d+u ),ir u+1 , ...,ir t ,p. 
From this structure, construct the following invariant tensor, 

m 1 ...m t n ._ m 1 ...m d m' m d+1 ...m u m" n m u+1 ...m t m (a nr\ 

C •— C(rf) C( u ) C(d+u)m'm" C(t)n j l 4 - Z0 J 

associated with the representations 7Ti, ...,n t ,p. To obtain a non-trivial result in the end, 
we need all the tensors to be non-zero. The existence of such tensors is equivalent to the 
following two conditions on the data (a-c): 

d) The representations p^) and //(„) emerge respectively, in ni®...®7i d and 77^+1 <8>...<8>7r u ; 
and, 

e) the representation p(d+u) emerges both in p(d) <8> P(u) and n u+ i <g> ... <g> ir t <S> p- 
Finally, introduce an extended spin network (7,7?, p,c,M) such that 

7? = (7Ti,...,7r t , ...,7rjv), p= (p,p 2 ,...,pv), c= (c,c 2 ,...,c y ), (4.26) 

the remaining entries being arbitrary. Then, the corresponding state M s1 q is an eigen- 
vector of the operators (J s ^l) 2 , (J^l) 2 , and {jf^) 2 with the eigenvalues + 1), 
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_|_ i) an d j(d+u)^j(d+u) _j_ -g^ respectively, where the half integers j^- d \j^ u ' and 
j(d+u) correS p 0nc l to the representations P(d),P( u ) and P(d+u)- Hence, this N is also 



an eigenvector of with the eigenvalue (|4.4j , [4.5|) . It is obvious, that for any triple of 
representations and P(d+u) satisfying the constraint ( |4.5| ) there exists an extended 

spin network (|4.26|) . 

This construction provides all eigenvectors of Ag )t ,. The key reason behind this com- 
pleteness is that, given any choice of 7i"i, iid, 7r u , n t and p as above, the invariant 
tensors which can be written in the form Q4.25| ) with any p J t d )i and fAu+u) span the 
entire space C v of invariant tensors at v compatible with that data. Since the defining 
formula for a spin network function (|4.19|) is linear with respect to every component of 
c, given any spin network (7, tt, p, c, M) it suffices to decompose the component C\ of c at 
v 1 = v into invariant tensors of the form ( 4.25Q in any manner to obtain a decomposition 
of the corresponding spin network function into a linear combination of extended spin net- 
work functions given by ( |4.25| , [4.26 ). The desired result now follows from the orthogonal 
decomposition of TL° in to the extended spin network subspaces. 

Let us now turn to the operator A$. A basis of eigen vectors can be obtained in the 
following way Since the area operator can be expressed in terms of and commutes with 
(J^]) 2 , (j's'l) 2 , and (J^v^) 2 at any point v in S, we can simultaneously diagonalize all 
these operators. Because for every graph the area operator preserves the subspace of spin- 
network states associated with that graph and for two different graphs the spin network 
spaces spaces are orthogonal, it is enough to look for eigen vectors for an arbitrary graph 
7. Given a graph 7, labelings 7?, p and M as in section 4.2, at every vertex v contained 
in the surface S choose a basis in the space C v consisting of invariant tensors of the form 
( |4.25| ). The set of the spin network functions (4.19) constructed by varying 7, 7?, p, M and 
picking at each vertex v an element of the basis in C v constitutes a basis in 7i°. (If we 
restrict the labelings to p consisting only of the trivial representations, then the resulting 
set of spin network states provide a basis for the space 7i° of gauge invariant functions.) 
Each of such states is automatically an eigenvector of A$ with eigenvalue (|4.6|) . 

We conclude the first part of this sub-section with a simple example of an eigen vector 



of the area operator with eigenvalue as, where as is any real number satisfying ( |4.6| . [4.6|) . 

• Example. Suppose {jffl,jj?\j£ +u '), « = 1, W, is a finite set of triples of half 
integers which for every a satisfy ( |4.5| )). Rather than repeating the construction a)- 
e) above step by step, we will specify only the simplest of the resulting (extended) 
spin-networks. In S choose W distinct points v a , a = 1, ...,W. To every point 
v a assign two finite analytic curves e^ a and e UjQ starting at v a , not intersecting 
S otherwise, and going in opposite directions of S. For a graph 7 take the graph 
{ed,i,e U) i, ■■■,ed,w,£u,w}, the vertices being the intersection points v a and the ends 
of the edges e^a and e u>a (the curves being chosen such that the points v a are the 
only intersections). Label each edge ed, a with the irreducible representation 7rd, a 
corresponding to a given and every edge e UjQ with the irreducible representation 
7r U)Q defined by That defines a labeling 7? of 7. (The absence of edges e t>Q is 
equivalent to introducing these edges in any manner and assigning to them the trivial 
representations.) To define a labeling p at the vertices v a , assign to every vertex 
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v a a representation p a defined by a given j^ +u \ Next, to each vertex v a assign an 
invariant tensor c™ dm " m associated to the triple of representations (7T</ )Ct , Ti u ,a, Pa) 
introduced above. The construction of a spin-network is completed by: (z) labeling 
that end point of each e^ a and, respectively, of e UjCe which is not contained in S, with 
the representation p^ a '■= ^d,a and, respectively, p u<a := 7T U}a ] (ii) labeling of these 
ends of the edges with the unique invariants corresponding to the representations 
H(d\ >a , P(d), a or ) respectively, to P( u ), a , P(u),a\ (^) defining a labeling M of vertices 
which can be chosen arbitrarily, provided at a vertex v a the associated vector M a 
belongs to the representation P(d+u),a and at an endpoint of either of the edges ed/ u ,a 
the associated M d / U , a belongs to ptd/u),a- 

As we noted in section 2, the Hilbert space 7i° is the quantum analog of the full phase 
space. Now, in the classical theory, the imposition of the Gauss constraint on the phase 
space does not restrict the allowed values of the functional As of (|2.2|) . It is therefore of 
interest to see if this feature persists in the quantum theory: Is the spectrum of A$ on 
the full 7i° the same as that on its gauge invariant sub-space 7i°? As was indicated in 
Sec. 4.1, the answer is in the affirmative only if the surface is open. If S is closed, there 
are restrictions on the spectrum which depend on topological properties of S embedded 
in E. The second part of this section is devoted to this issue. As indicated in Sec. 4.1, 
we need to consider three separate cases. 

Case i): dS ^ (and dS cEj. 

We will modify the spin-network of the above Example in such a way as to obtain 
a gauge-invariant eigenstate without changing the eigen value of the area operator. Let 
7 and the labeling jf be the ones defined in the Example. To each vertex v a assign one 
more edge e t>a beginning in v a and contained in S. Label it by the representation n t >a 
corresponding to a given j^ +uS> at that point. The labeling p is now taken to be trivial. 
To every point v a assign, as in the Example, an invariant tensor c a associated now to the 
representations {^d,a, ftu,a, ^t,a)- Every extension of this data to a spin-network will define 
a spin-network state which is gauge invariant at each of the points v a . Now, we need to 
define a closed spin-network which contains all the edges e^a, e UjCC , e tt a and provides an 
extension for the labelings already introduced. For this, we use a key property of the 
area operator associated to a surface with boundary: vertices which lie on dS do not 
contribute to the action of the operator. Therefore, we can simply extend every edge e tya 
within S to the boundary of S. Denote the intersection point with dS by Vt >a - Next, for 
every a we extend (in a piecewise analytic way) the edges ed, a and e UjQ such that they end 
at v ti a- The extended edges form a graph 7' = {e'^i, e' Uj i, e'^i, e'd,w, z'u,w, e't,w}- Let 
us label each primed edge by the irreducible representation assigned before to the edge 
it is an extension of. This defines a labeling ir' of 7'. Finally, assign to each new vertex 
v ti0l the non-zero invariant tensor c' t ,am u ,m d ,m t (which is unique up to rescaling) associated 
to the triplet of representations (7id,a, ^u,a, ^t,a)- This completes the construction of a 
gauge invariant extension of a spin-network state constructed in the Example. Thus, for 
an open surface, the spectrum of the area operator A s on H° is the same as that on 7i°. 
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Case ii): dS = and S splits £ in to two open sets. 

In this case we can not repeat the above construction: Since S has no boundary, if 
additional vertices are needed to close the open spin-network, they must now lie in S and 
can make unwanted contributions to the action of the area operator. Consequently, there 
are further restrictions on the possible eigenvalues of the operators {j'sl) 2 , {jj^l) 2 and 
(Jsv ) 2 - To see this explicitly, consider arbitrary spin-network state (7,7?, c) given by 
the construction a)-e) of Sec 4.3. Let {vi, ...,vw} be a set of the vertices of 7 contained 
in the surface S. Graph 7 can be split into three graphs: 7 t which is contained in S, 
7 n which is contained in a one side of S in E and 7^ contained in the other side of S 
in S. The only intersection between the two parts is the set {v%, ...,vw} of vertices of 
7 which are contained in S. Let 7 r be one of the parts of 7 (i.e. r = d or r = u or 
r = t). According to the construction a)-e), the labelings 7? and c define naturally on 7 r 
an extended spin-network. The labeling of the edges of 7 r by irreducible representations 
is defined just by the restriction of 7? to j r . The labeling of the vertices by irreducible 
representations and invariant tensors is defined in the following way. For the vertices of 
7 r which are not contained in S, the labelings are taken to be again the restriction of p 
(which are all trivial) and c. To a vertex v a contained in S we assign the representation 
corresponding to a given j riCt and the invariant tensor c r defined in b) (for r = d,u) and 
c) (for r = t) of the construction a) - e). Finally, we complete it by arbitrary nonzero 
labeling M of the vertices with vectors in appropriate representations. The construction 
a) - e) guarantees that a resulting extended spin- network state is not zero. Now, for an 
extended spin-network (7', jf, p ', cf, M') we have the following "fermion conservation law" : 



EVW =N (4.27) 

V 

for some integer N, where v runs through the vertices of a graph 7' and each is 
an half-integer corresponding to an representation assigned to v by p . In our case we 
therefore obtain the restriction: 

J2jr,a =N r (4.28) 

a 

for r = d,u,d + u which gives the conditions (|4.7|) listed in Sec 4.1. (In fact either two of 
the above conditions imply the third one). 

The conditions ( |4.28| ) are also sufficient for an eigen vector to exist. Suppose we are 
given a set of half integers as in the Example above, which satisfy the restriction ( |4.28|) . 
A statement "converse to the fermion conservation law" is that for any set {t>i, ...,vw} of 
points in S and any assignment v a h-> j a where j a are non-negative half integer satisfying 
( [L28j ), there exists an extended spin-network (7', 7?', p, C?, M') such that every v a is its 



index, j a corresponds to the representation assigned to v a by p, and for every vertex 
v 7^ v a , a = 1,...,W, of 7', the representation assigned by p is trivial. From extended 
spin-networks provided by the above statement it is easy to construct an eigen vector of 
the corresponding eigen values. 
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Case Hi): dS = but S does not split S. 

The only difference between this case and the previous one is that now a graph 7 
representing an eigen vector is cut by S into two components: j t contained in S and 
■jd+u which corresponds to the rest of 7. Since jd+u can be now connected by the same 
arguments as above, we prove that a necessary and sufficient condition for an eigen vector 
to exists is (|4.28|) imposed only on the half integers +u >. 



V. DISCUSSION 

In section 1, we began by formulating what we mean by quantization of geometry: 
Are there geometrical observables which assume continuous values on the classical phase 
space but whose quantum analogs have discrete spectra? In the last two sections, we 
answered this question in the affirmative in the case of area operators. In the next paper 
in this series we will show that the same is true of other ('3-dimensional' ) operators. The 
discreteness came about because, at the microscopic level, geometry has a distributional 
character with 1-dimensional excitations. This is the case even in semi-classical states 
which approximate classical geometries macroscopically [0,0 • 



We will conclude this paper by examining our results on the area operators from 
various angles. 

1. Inputs: The picture of quantum geometry that has emerged here is strikingly 
different from the one in perturbative, Fock quantization. Let us begin by recalling the 
essential ingredients that led us to the new picture. 

This task is made simpler by the fact that the new functional calculus provides the 
degree of control necessary to distill the key assumptions. There are only two essential 
inputs. The first assumption is that the Wilson loop variables, T a = TrV — exp J a A, 
should serve as the configuration variables of the theory, i.e., that the Hilbert space of 
(kinematic) quantum states should carry a representation of the C*-algebra generated 
by the Wilson loop functionals on the classical configuration space A/Q. The second 
assumption singles out the measure p,°. In essence, if we assume that Ef be represented by 
— ih5/5A l a , the 'reality conditions' lead us to the measure jl° QTUJ . Both these assumptions 
seem natural from a mathematical physics perspective. However, a deeper understanding 
of their physical meaning is still needed for a better understanding of the overall situation 

B 

Compactness of SU(2) plays a key role in all our considerations. Let us therefore 



briefly recall how this group arose. As explained in [0,pl|, one can begin with the ADM 



8 In particular, in the standard spin-2 Fock representation, one uses quite a different algebra 
of configuration variables and uses the flat background metric to represent it. It then turns 
out that the Wilson loops are not represented by well-defined operators; our first assumption is 
violated. One can argue that in a fully non-perturbative context, one can not mimic the Fock 
space strategy. Further work is needed, however, to make this argument water-tight. 
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phase-space in the triad formulation, i.e., with the fields {Ef, K l a ) on S as the canonical 
variables, and then make a canonical transformation to a new pair (A l a := (r* + K l a ), Ef), 
where K l a is the extrinsic curvature and T l a , the spin-connection of Ef. Then, A l a is an 
SU(2) connection, the configuration variable with which we began our discussion in sec- 
tion 2. It is true that, in the Lorentzian signature, it is not straightforward to express 
the Hamiltonian constraint in these variables; one has to introduce an additional step, 
e.g., a generalized Wick transform [IS]. However, this point is not directly relevant in 
the discussion of geometric operators which arise at the kinematical level. (See, however, 
below). Finally, we could have followed the well-known strategy of simplifying con- 
straints by using a complex connection 9A* := (T^ — in place of the real A\. The 
internal group would then have been complexified 577(2). However, for real (Lorentzian) 
general relativity, the kinematic states would then have been holomorphic functionals of 
A l a . To construct this representation rigorously, certain technical issues still need to be 
overcome. However, as argued in |ISJ, in broad terms, it is clear that the results will be 
equivalent to the ones obtained here with real connections. 

2. Kinematics versus Dynamics: As was emphasized in the main text, in the classical 
theory, geometrical observables are defined as functionals on the full phase space; these 
are kinematical quantities whose definitions are quite insensitive to the precise nature of 
dynamics, presence of matter fields, etc. Thus, in the connection dynamics description, 
all one needs is the presence of a canonically conjugate pair consisting of a connection and 
a (density weighted) triad. Therefore, one would expect the result on the area operator 
presented here to be quite robust. In particular, they should continue to hold if we bring 
in matter fields or extend the theory to super-gravity. 

There is, however, a subtle caveat: In field theory, one can not completely separate 
kinematics and dynamics. For instance, in Minkowskian field theories, the kinematic field 
algebra typically admits an infinite number of inequivalent representations and a given 
Hamiltonian may not be meaningful on a given representation. Therefore, whether the 
kinematical results obtained in any one representation actually hold in the physical theory 
depends on whether that representation supports the Hamiltonian of the model. In the 
present case, therefore, a key question is whether the quantum constraints of the theory 
can be imposed meaningfully on 7Y°|] Results to date indicate (but do not yet conclusively 
prove) that this is likely to be the case for general relativity. The general expectation is 
that this would be the case also for a class of theories such as super-gravity, which are 
'near' general relativity. The results obtained here would continue to be applicable for 
this class of theories. 

3. Dirac Observable: Note that As has been defined for any surface S. Therefore, 
these operators will not commute with constraints; they are not Dirac observables. To 
obtain a Dirac observable, one would have to specify S intrinsically, using, for example, 



Note that this issue arises in any representation once a sufficient degree of precision is reached. 
In geometrodynamics, this issue is not discussed simply because generally the discussion is rather 
formal. 
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matter fields. In view of the Hamiltonian constraint, the problem of providing an explicit 
specification is extremely difficult. However, this is true already in the classical theory. In 
spite of this, in practice we do manage to specify surfaces and furthermore compute their 
areas using the standard formula from Riemannian geometry which is quite insensitive 
to the details of how the surface was actually defined. Similarly, in the quantum theory, 
if we could specify a surface S intrinsically, we could compute the spectrum of As using 
results obtained in this paper. 

4. Comparison: Let us compare our methods and results with those available in 
the literature. Area operators were first examined in the loop representation. The first 
attempt |12j was largely exploratory. Thus, although the key ideas were recognized, the 
very simplest of loop states were considered and the simplest eigenvalues were looked 
at; there was no claim of completeness. In the present language, this corresponds to 
restricting oneself to bi-valent graphs. In this case, apart from an overall numerical factor 
(which does, however, have some conceptual significance) our results reduce to that of 
0- 

A more complete treatment, also in the framework of the loop representation, was 
given in |]i~3| . It may appear that our results are in contradiction with those in ]13j on 



two points. First, the final result there was that the spectrum of the area operator is 
given by £ 2 P J2 \l + 1) 5 where ji are half-integers, rather than by ( |4.6| ). However, the 
reason behind this discrepancy is rather simple: the possibility that some of the edges at 
any given vertex can be tangential to the surface was ignored in [13] . It follows from our 
remark at the end of section 4.2 that, given a surface S, if one restricts oneself to only to 
graphs in which none of the edges are tangential, our result reduces to that of ||13|| . Thus, 
the eigenvalues reported in [T3J do occur in our spectrum. It is just that the spectrum 
reported in |13| is incomplete. Second, it is suggested in [LpJ that, as a direct consequence 
of the diffeomorphism covariance of the theory, local operators corresponding to volume 
(and, by implication, area) elements would be necessarily ill-defined (which makes it 
necessary to by-pass the introduction of volume (and area) elements in the regularization 
procedure). This assertion appears to contradict our finding that the area element yj~gs 
is a well-defined operator-valued distribution which can be used to construct the total 
area operator A$ in the obvious fashion. We understand [29], however, that the intention 
of the remark in [|T3J was only to emphasize that the volume (and area) elements are 
'genuine' operator-valued distributions; thus there is no real contradiction. 



The difference in the methodology is perhaps deeper. First, as far as we can tell, in [|L3 



only states corresponding to trivalent graphs are considered in actual calculations. Thus, 
even the final expression (Eq (48) in [|H|) of the area operator after the removal of the 
regulator is given only on tri-valent graphs, similarly, their observation that every spin 
network is an eigenvector of the area operator holds only in the tri-valent case. Second, for 
the limiting procedure which removes the regulator to be well-defined, there is an implicit 
assumption on the continuity properties of loop states (spelled out in detail in |30|). A 
careful examination shows that this assumption is not satisfied by the states of interest 
and hence an alternative limiting procedure, analogous to that discussed in section 3.1, is 
needed. Work is now progress to fill this gap fl29| . Finally, not only is the level of precision 
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achieved in the present paper significantly higher but the approach adopted is also more 
systematic. In particular, in contrast to ||13|| , in the present approach, the Hilbert space 
structure is known prior to the introduction of operators. Hence, we can be confident 
that we did not just omit the continuous part of the spectrum by excising by fiat the 
corresponding sub-space of the Hilbert space. 

Finally, the main steps in the derivation presented in this paper were sketched in the 
Appendix D of [ID|. The present discussion is more detailed and complete. 

5. Manifold versus Geometry: In this paper, we began with an orientable, analytic, 
3-manifold S and this structure survives in the final description. As noted in footnote 
1, we believe that the assumption of analyticity can be weakened without changing the 
qualitative results. Nonetheless, a smoothness structure of the underlying manifold will 
persist. What is quantized is 'geometry' and not smoothness. Now, in 2+1 dimensions, 
using the loop representation one can recast the final description in a purely combinatorial 
fashion (at least in the so-called 'time-like sector' of the theory). In this description, at 
a fundamental level, one can avoid all references to the underlying manifold and work 
with certain abstract groups which, later on, turn out to be the homotopy groups of the 
'reconstructed/derived' 2-manifold (see, e.g., section 3 in |31]|). One might imagine that, 
if and when our understanding of knot theory becomes sufficiently mature, one would also 
be able to get rid of the underlying manifold in the 3+1 theory and introduce it later as 
a secondary /derived concept. At present, however, we are quite far from achieving this. 

In the context of geometry, however, a detailed combinatorial picture is emerging. 
Geometrical quantities are being computed by counting; integrals for areas and volumes 
are being reduced to genuine sums. (However, the sums are not the 'obvious' ones, 
often used in approaches that begin by postulating underlying discrete structures. In the 
computation of area, for example, one does not just count the number of intersections; 
there are precise and rather intricate algebraic factors that depend on the representations 
of SU(2) associated with the edges at each intersection.) It is striking to note that, in the 
same address |3[] in which Riemann first raised the possibility that geometry of space may 
be a physical entity, he also introduced ideas on discrete geometry. The current program 
comes surprisingly close to providing us with a concrete realization of these ideas. 
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